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Abstract. We investigate the geodesics in the entire class of nonexpanding impulsive grav¬ 
itational waves propagating in an (anti-)de Sitter universe using the distributional form of 
the metric. Employing a 5-dimensional embedding formalism and a general regularisation 
technique we prove existence and uniqueness of geodesics crossing the wave impulse leading 
to a completeness result. We also derive the explicit form of the geodesics thereby confirming 
previous results derived in a heuristic approach. 
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1. Introduction 


Impulsive pp-waves have now been studied for several decades and have become textbook 
examples of exact radiative spacetimes modelling short but intense bursts of gravitational 
radiation propagating in a Minkowski background (see e.g. [ISl Sec. 20.2] and the references 
therein). Such geometries have been introduced by Roger Penrose using his ‘scissors and 
paste method’ (see e.g. m) leading to the distributional Brinkmann form of the metric 

(1.1) ds^ = dx^ -h dy^ — 2dudv -|- f{x, y)6{u)du‘^, 

i.e., impulsive limits of sandwich pp-waves [5]. Alternatively, almost at the same time, Aichel- 
burg and Sexl in [T] have obtained a special impulsive pp-wave as the ultrarelativistic limit 
of the Schwarzschild geometry and several authors have applied the same approach to other 
solutions of the Kerr-Newman family (see e.g. [H Ch. 4] and [271 Sec. 3.5.1] for an overview). 
This procedure of boosting static sources to the speed of light was later generalised to the case 
of a non-vanishing cosmological constant A in the pioneering work m by Hotta and Tanaka 
(see also |28l [29] ) which lead to an increased interest in nonexpanding impulsive waves in cos¬ 
mological de Sitter and anti-de Sitter backgrounds. The Penrose ‘scissors and paste method’ 
for non-vanishing A was described in |38l [30] while impulsive limits in the Kundt class were 
considered in |26| and elsewhere, see m Sec. 20.3] for an overview. 


Generally, nonexpanding impulsive waves in all backgrounds of constant curvature can be 
described by a continuous as well as by a distributional form of the metric tensor. To give a 
brief discussion of these we start with the conformally flat form of Minkowski (A = 0) and 
(anti-)de Sitter (A 0) background spacetimes 


( 1 . 2 ) 


2 dry dy — 2 dW dV 
[l+lA(ryfy-7/V)]2- 
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As in m. here U, V are the usual null and r], rj the usual complex spatial coordinates. Now 
for lA > Q WQ perform the transformation 


(1.3) 


U : 


■-U, V = V + H + UH^zH^z , 


where H{Z, Z) is an arbitrary real-valued function. Combining this with the background line 


element (1.2) for [/ < 0 in which we set lA = U, V = V, rj = Z we obtain the continuous form 


of the metric 


(1.4) 


2 2\dZ + U+{Hzz^Z + Hzz<^Z)\^ -‘^^UdV 
ds = — 


[1 + I.K{ZZ -UV + U+G)\^ 


where]] G{Z, Z) = ZH^z + ZHz - H and U+ = U+{U) = 0 for [/ < 0 and U+{U) = U for 
?7 > 0. This ‘kink-function’ 17+ is Lipschitz continuous, hence the spacetime (apart from 
possible poles of H, which indeed occur in physically realistic models, see e.g. mm, and 
Section below) is locally Lipschitz. Recall that a locally Lipschitz metric possesses a 
locally bounded connection and hence a distributional curvature, which, however, in general 
is unbounded. In fact the discontinuity in the derivatives of the metric introduces impulsive 
components in the Weyl and curvature tensors (see m), and the metric (|1.4|) explicitly 


describes impulsive waves in de Sitter, anti-de Sitter or Minkowski backgrounds. 


For A = 0, (1.4) reduces to the classic Rosen form of impulsive pp-waves (cf. |16( Sec. 17.5]). In 
this special case the geodesic equation has been rigorously solved in |23] using Caratheodory’s 
solution concept (see e.g. HI Ch. 1]), which allows to deal with the locally bounded but 
discontinuous right hand side of the equation. The geodesics thereby obtained coincide with 


the limits of the geodesics for the distributional form (1.1) calculated in m which have been 
previously derived heuristically (e.g. in [TTl [2l [35]L 

To deal, however, with the geodesic equation for the full class of nonexpanding impulsive waves 
for arbitrary A, that is the complete metric 0. the more sophisticated solution concept 
of Filippov m Ch. 2]) has been applied recently in |32]. Building on a general result 
for all locally Lipschitz spacetimes (|lQ]), existence, uniqueness and global C^-regularity of 
the geodesics has been established. This, in particular, justihes the C^-matching procedure 
which has been used before to explicitly derive the geodesics in this and similar situations 

([saSHEniEIllMllM!). 

On the other hand, the distributional form of the impulsive metric arises by writing the 


transformation relating (1.2) and (1.4) in a combined way for all U using the Heaviside 
function 0 as 


(1.5) 


U = u, V = V + eH + U+H^zHz, r] = Z + U. 


+ H^z- 


Of course, (1.5) is discontinuous in the coordinate V and merely Lipschitz continuous in rj 


across {lA = 0} but applying it formally to (1.4) we arrive at 

2 dr/ dr? - 2 dZ7 dV + 277(r/, ry) 6{U) 


( 1 . 6 ) 


ds^ = 


T + iA(r/r7-L/V)]2 


which has the striking advantage of coinciding with the background metric dsg (1.2) off the 
impulse located at lA = 0. This, however, comes at the price of introducing a distributional 
coefficient in the metric which leads us out of the Geroch-Traschen class m ) of metrics (of 


^This choice of sign of G is in accordance with |31| . which is our main point of reference, but different from 
more recent papers, e.g. m- 
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which guarantees existence and stability of the curvature in distri- 


regularity n 

butions (see also [Mil]). However, due to its simple geometrical structure the metric ( |1.6| ) 
nevertheless allows to calculate the curvature as a distribution, again leading to impulsive 
components in the Weyl and curvature tensors i |30j 1. Also, in Minkowski background the 


metric (1.6) reduces to the Brinkmann form of impulsive pp-waves (1.1). 


Clearly, a mathematically sound treatment of the transformation (1.5) is a delicate matter. 
In the special case of pp-waves this has been achieved in 


using nonlinear distributional 
geometry ([m Ch. 3]) based on algebras of generalised functions ([H])- More precisely, the 
‘discontinuous coordinate change’ was shown to be the distributional limit of a ‘generalised 
diffeomorphism’, a concept further studied in mm- The key to these results was provided 


by a nonlinear distributional analysis of the geodesics in the metric (1.1), and, in particu¬ 


lar, an existence, uniqueness and completeness result for the geodesic equation in nonlinear 
generalised functions mm)- 

This suggests that a first step towards the long-term goal of understanding the transforma¬ 


tion (1.5) for A 7 ^ 0 is to reach a mathematically sound understanding of the geodesics in 
the distributional metric (1.6). The geodesic equation for (1.6), however, displays a very 
singular behaviour including terms proportional to the square of the Dirac-J (cf. |38]1. For 
this reason the authors of |31j have employed a five-dimensional formalism ( |29L [30] ) of em¬ 
bedding (anti-)de Sitter space into a 5-dimensional pp-wave spacetime (see Section below). 
In this approach the geodesic equation takes a form that is distributionally accessible at all, 
however, not mathematically rigorously. In the absence of a valid solution concept for this 
nonlinear distributional equations a natural ansatz was used to derive the geodesics and to 
study them in detail in m Sec. 4-5]. Nevertheless, a desirable nonlinear distributional anal¬ 
ysis of the geodesic equation in the A ^ 0-cases, which will eventually lead to a mathematical 


understanding of the transformation (1.5), has been missing to date. 

In this paper we provide such an analysis. Thereby, we follow [31] in employing the five¬ 
dimensional formalism. We will, however, not use any theory of nonlinear distributions leaving 
a detailed study of nonexpanding impulsive waves in (anti-)de Sitter universe as (nonlinear) 
distributional geometries to a subsequent paper. Instead, we will employ a regularisation 


approach and view (1.6) as a spacetime with a short but hnitely extended impulse (i.e., a 
generic sandwich wave with support in a regularisation strip which we will also call the ‘wave 
zone’) and employ an analysis in the spirit of |36] where impulsive limits in a class of pp-wave 
type spacetimes with a more general wave surface but vanishing A ( (6] HS] [T] E] ) have been 
considered. 

We will detail this regularisation approach in the next section after introducing the 5-dimen¬ 


sional formalism.In particular, we will replace the Dirac-5 in the metric (2.1) below by a very 
general regularisation 5^, thereby ensuring that our results are regularisation independent 
(within the class of so called a model delta nets). Then in Section]^ we will employ a hxed point 
argument to show that the regularised equations have unique smooth solutions which cross 
the regularisation strip. This will lead to our main result on completeness of nonexpanding 
impulsive gravitational waves in a cosmological background. The technical proofs allowing 
for the application of the fixed point theorem are deferred to Appendix In Section we 
study boundedness properties of the regularised geodesics which are essential when dealing 
with their limits in Section There we show that the solutions of the regularised geodesic 
equation converge, as the regularisation parameter goes to zero, to geodesics of the background 
(anti-)de Sitter spacetime which have to be matched appropriately across the impulse and 
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have been derived previously in m- Overly technical calculations are collected in Appendix 


2 . The geodesic equation for A / 0 

In this section we detail our regularisation approach and derive the respective geodesic equa¬ 
tions. To begin with, however, we recall the 5-dimensional formalism of |29l [30] for the full 
class of nonexpanding impulsive waves with non-vanishing A. To this end we start with the 


5-dimensional impulsive pp-wave spacetime M with metric (extending (1.1)) 


( 2 . 1 ) 


ds^ = dZ| + dZ| + adZl - 2dUdV + H{Z2, Z^, Zi)5{U)dU‘^ 


and consider the four-dimensional (anti-)de Sitter hyperboloid {M,g) given by the constraint 

( 2 . 2 ) + aZj - 2UV = aa^, 

where a := (crA), cr ;= sign(A) = ±1, and U = + Zi ), V = :^(-^0 ~ -^i) are null- 

coordinatesj^ Here (Zq, ..., Z 4 ) are global Cartesian coordinates of M^. The impulse is located 
on the null hypersurface {[/ = 0 } given by 

( 2 . 3 ) Z2^ Z^^ oZ/^ = 00 ? , 


which is a nonexpanding 2-sphere in the de Sitter universe (A > 0) and a hyperboloidal 2 - 
surface in the anti-de Sitter universe (A < 0), respectively. Various 4-dimensional coordinate 
parametrizations of these spacetimes have been considered e.g. in |28j . Physically the space¬ 


time ( 2 . 1 ), ( 2 . 2 ) describes impulsive gravitational waves as well as impulses of null matter. 
Purely gravitational waves occur in case the vacuum field equations are satished. It was 
demonstrated in [29l [30] that such solutions can be explicitly written as 

(2.4) 


H{z,4>) = boQl{z) + J2m=l ^rn {z) COs[m{(t) - ^m)] 


where z = ^ 4 / 0 , tancj) = Z 3 /Z 2 and Q^{z) are associated Legendre functions of the sec¬ 
ond kind generated by the relation Q'{^{z) = (—cr)”*|l — -^^Qi{z). The first term for 


m = 0, i.e., Qi{z) = § log 


1 - 1-2 


1—Z 


— 1, represents the simplest axisymmetric Hotta-Tanaka so¬ 


lution m)- The components with m > 1 describe nonexpanding impulsive gravitational 
waves in (anti-)de Sitter universe generated by null point sources with an m-pole structure, 
localized on the wave-front at the singularities 2 ; = ± 1 . See [aolEHlEZI for more details. 


Now, the geodesics 7 of M with tangent T are characterized by the condition that their M- 
acceleration A = V^T is everywhere normal to M. Denoting by N the normal vector to M 
in M with g{N, N) = cr, we hence obtain 


(2.5) 


VtT = -ag{T,VTN)N. 


^These coordinates are different from those used in the metric (1.4|. Since in this paper we will not use the 
continuous form (1.41 we simplify the notation by not distinguishing them by a bar (as we did in [32]). 
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Using this identity and the constraint (2.2) the explicit form of the geodesic equation was 
given in [311 eq. (28)] as 

U = -^-AUe, 

V -\h 5'U^ -5P^H^p8ZqU = -\KV{e+\G5U^) , 

= -^KZi{e + ]^G5U‘^), 

= -^AZ^{e+^G6U^) , 


( 2 . 6 ) 


where 

(2.7) 


G :=5P'i ZpH^q- H, and e := ff(r,T) = ±1, 0 


denotes the normalisation constant for spacelike (e = 1), timelike (e = —1) and null (e = 0) 
geodesics respectively. Observe, that' denotes the derivative with respect to an affine param¬ 
eter t which we have suppressed in the equations. Here and in the following we also adopt 
the convention that Greek indices a, f3 take all values 0,1,2,3,4 while the indices p,q,r are 
restricted to the values 2, 3,4 and j run from 2 to 3 only. 

Obviously, these equations reduce to the geodesic equations of the (anti-)de Sitter background 
off the impulse located at {U = 0}. Also observe that the first equation decouples from the 
rest of the system and can be easily integrated. Consequently U can be used as a parameter 


of the remaining equations, a fact which is essential for the analysis of the system (2.6) 


presented in ini Sec. 4]. In fact, there the geodesics of the (anti-)de Sitter background in 
front and behind the wave impulsive are matched using a natural ansatz for solutions in the 
entire spacetime. However, this procedure has to be viewed as being only heuristic since the 
solution’s Zp-components are (assumed to be) continuous but not C^, while the U-component 
is (assumed to be) even discontinuous. Consequently the solutions cannot be plugged back 
into the equations due to the occurrence of undefined products of distributions, and so the 
question arises in which sense they actually solve the equations, see the discussion at the end 
of Sec. 4 of [31j . The situation is similar to the one encountered for impulsive pp-waves with 
A = 0 and we refer to the discussion in [391 Sec. H] as well as to the general discussion in 

m. 


To resolve this open problem we now employ a regularisation approach and detail the setting 


we are working with: To begin with we replace the Dirac-5 in the line element (2.1) by a fairly 


general class of smooth approximations called model delta nets. Chose an arbitrary smooth 
function /? on M with unit integral and its support contained in [—1,1]. Then for 0 < e < 1 
set 

( 2 . 8 ) 5e{x) ■.= ^^p{^ . 

We now consider for fixed e G (0,1] the five-dimensional sandwich wave 

(2.9) ds^ = dZl + dZl + adZl - 2dUdV + 77(^2, Z3, Zi)8^{U)dU^ , 


and define the spacetime of our interest as {M,g^) given by the constraint (2.2), i.e., 
(2.10) F{U, U, Z2, Z3, Z 4 ) := -2UV + - aa^ = 0 . 
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Observe that while the differential dF = 2(—1/, —[/, Z2, Z3, (TZ4) is independent of e, the 
normal vector iV" := g^^dFp depends on e. Indeed we choose to work with the non-normalised 
normal vector to M given by 

(2.11) Ne = {U,V + HU6eiU),Zp) with g,{N,,N,) = aa^ -U^H5e{U). 

The non-zero Christoffel symbols of ( |2.9| ) are given by 

( 2 . 12 ) T^,uu = -\h6',{U) , = -\h,p5,{U) , 

(2.13) Tluu = -\hMU) , , 

and the geodesics of {M,g^) are now characterised by 

(2.14) = 

where (suppressing the parameter) we write yg = (14,14, Zp^) for the geodesics with tangent 
Tg = {U^,Vs, Zp^) and V denotes the Levi-Civita connection of (2.9). By a straightforward 
calculation we now obtain 

1 


g,iT„ V^.iVg) = e + ^U^ Ge{Ue, Zp,) - U, (H{Zp,) 6,{U,) U, 


(2.15) 

where we have used the abbreviations 

(2.16) G,{U,Zr):=6P'iH^p{Zr)5e{U)Zq + H{Zr)5',{U)U , and e := ^.(Tg, Tg) = ±1,0 . 

Observe that since the g'g-norm of the tangent vector Tg is constant along the geodesic 7g, we 
have chosen it also to be constant in e, which means that we have fixed the normalisation 
independently of e. Finally we obtain the following explicit form of the geodesic equations 


ii. = -(e + ^ f/g" Ge - Ue {H6e U, 


1/g 


Ve--H6, 1/2 _ spm,p 6,ZlUe = -(e + - 1/| Gg - 1/g {H 6, U, 


(2.17) 


5eUe = - (e + - 1/2 Gg - 17g (77 4 U, 


Z^e - ^ 77,4 417' = -(e + \ulG,-Ue{H5e U, 


(Ta^ 

- UIH5, 

Hg + 7741 /e 

aa? 

- UIH5, 


Zis 

aa? 

- UlH5e 


Zie 

aa? 

- UlH6e 


where we again have suppressed the parameter t as well as the dependencies on the variables. 
However, note that always 

4 = 4(t7g(t)), 5( = (5((l/g(t)), 

(2.18) G, = G,iU,{t),Zp,{t)), 77 = 77(Zpg(t)), and Hp = Hp{Zg,it)). 

The right hand sides of these equations are considerably more complicated than their ‘dis¬ 
tributional counterparts’ in (2.6), the reason being that in the regularised equations the 
distributional identities (5(17)1/ = 0 and 4(17)1/ = —(5(1/) do not apply. Indeed the lack of the 
first one leads to the more complicated form of the normal vector (see (2.11)) and is reflected 
in the second summand in the denominator of (2.17). On the other hand the lack of the 
second one is responsible for the different form of the terms involving Gg as compared to the 
ones involving G (cf. (2.7)) in (2.6), to which they reduce in the limit e —0. Finally the 
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terms proportional (if (JeC/g)’which occur in all four equations vanish in the limit e —)■ 0 again 
due to the first identity. The same holds true for the term proportional to contained in 

the hrst equation. In this sense the equations (2.17) converge weakly to (2.6) as e —)■ 0. 

The more complicated form of the system (2.17), in particular, results in the fact that the 
t4-equation does not decouple from the rest of the system and consequently U cannot be used 
as a parameter along the geodesics. This issue greatly complicates our analysis. However, 
the 14 -equation still is linear and decoupled, hence can be simply integrated once the rest of 
the system is solved. 


3. Existence, uniqueness and completeness of geodesics 


In this section we prove an existence and uniqueness result for solutions of the initial value 
problem for (2.17) that additionally guarantees that the geodesics that enter the sandwich 
region U G [—e, e] at one side exist long enough to leave it on the other side. This allows us 
to obtain global solutions of the geodesic equations, since outside the strip {—e < ?7 < s} the 
spacetime coincides with the background (anti-)de Sitter universe. 

Observe that for fixed e the equations are smooth and hence a local solution is guaranteed to 
exist. However, the time of existence might depend on e and in principle could even shrink 
to zero as e —)> 0. So the main objective here is to provide a result which guarantees that the 
time of existence is independent of e, and large enough such that the solutions pass through 
the regularisation sandwich (at least for all small e). To this end we employ a fixed point 
argument in the spirit of [36l Appendix A] based on Weissinger’s fixed point theorem f|33]l. 
However, the signihcant increase in the complexity of the equations forces the use of new ideas 
to derive the required estimates. In particular, since it is not possible to use the [/-coordinate 
as a parameter along the geodesics, the ‘singular terms’ such as 5^ are composed with the 
[/-coordinate of the solution see ( |2.18[ ). We have separated the technical proofs preparing 
the grounds for the application of the fixed point theorem from the main line of arguments 
and have deferred them to Appendix [A| 

Let us start by giving the general setup. Consider any geodesic 


(3.1) ^ = {U,V,Zp) 

on the background (anti-)de Sitter universe without impulsive wave but reaching [/ = 0. All 
other geodesics are not of interest to the present analysis and will be dealt with separately. We 
choose an affine parameter t in such a way that U{t = 0) = 0 and assume 7 to be normalised 
by e = ±1,0. Such geodesics can explicitly be written as 

(3.2) U = t, U = ail^ sinh{t / a), U = aU^ sm{t / a), 

in the cases ae = 0, ae < 0, and ere > 0, respectively, see lai eq. (33)]. Recall that the case 
fie = 0 corresponds to null geodesics in both de Sitter and anti-de Sitter space, while the 
case ae < 0 corresponds to timelike geodesics in de Sitter as well as to spacelike geodesics 
in anti-de Sitter spacetime, and hnally ae > 0 describes spacelike geodesics in de Sitter and 
timelike geodesics in anti-de Sitter space. Without loss of generality we assume the constant 
[7° to be positiv^so that in all three cases U is increasing (at least for t G [—a 7 r/ 2 , o 7 r/ 2 ]). 
It is thus most convenient to prescribe initial data at t = 0, that is we set 

(3.3) ^(i = o)=:(0,H°,zO), 7(t = 0) =: ([/° > 0, l/°, 4°), 


^The time reversed case with f/° < 0 can be treated in complete analogy. 
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where the constants satisfy the constraints 

(3.4) {Zlf + {Zlf + a{Zlf-2U^V^ = aa^, Z^Z^ +Z^Z^ + aZlzl-V^U^-U°V^ = 0 , 

which are simply consequences of the fact that we are dealing with 7 on the (anti-)de Sitter 
manifold, see (2.2). Note however that = 0, so the last term on the left hand side of either 
condition actually vanishes. In addition the normalisation condition 

(3.5) - 2t7°I>° + {Z^f + {Z^f + a{Zlf = e 
holds. 


Now we start to think of 7 as a geodesic in the impulsive wave spacetime (2.1), (2.2) ‘in front’ 
of the impulse that is for U < 0. Also, 7 is a geodesic in the regularised spacetime (2.9), 
(2.10) ‘in front’ of the sandwich wave, that is for U < —e. We will call it ‘seed geodesic’ and 
denote the affine parameter time when 7 enters this regularisation wave region by a^, 

(3.6) U{t = as) = —£. 

Observe that, by continuity of 7 , —?■ 0 from below as e —)• 0. More precisely, we have 

«£ = —e, ctg = —aArcsinh(e/al7°), and = —a arcsin(e/al7°), respectively for the three 
cases in (3.2), leading to + O(e^) in the latter cases and hence overall 

(3.7) - Ce < «£ < 0 , 
for some positive constant C. 


To investigate the geodesics in the regularised spacetime (2.9), (2.10), which is the main topic 
of this paper, we follow 7 up to the beginning of wave zone, i.e., up to t = a^, and then 
extend it (smoothly) to a geodesic 


(3.8) 


7e = (t/s, Fe, 


solving the regularised geodesic equations (2.17). This means 7 ^ at assumes the data (see 
Figure 1) 

(3.9) 7 e(«e) := 7(«.) =: (-£, Z %), 7 ,( 0 ,) := 7 ( 0 ,) =: . 

Observe that by smoothness of 7 the data ^eip-e) and 'fe{ois) converge to 7 ( 0 ) and 7 ( 0 ), 
respectively, as e —)• 0. In fact, by a mean value argument and we even have 

|(-e,v;°,zO)-(o,y° 




(3.10) 


\{UlV^,Z%)-{U^,V^,Z\ 


0 A 0 \ 


< 

sup || 7 (t)|U|a£| 

VI 

C£ 


ae<t<0 



< 

sup || 7 (t)||/i loel 

VI 

C£ 


Q;g:<t<0 




where h is any Riemannian background metric and C again is a generic constant. 

Based on Theorem A .6 proved in Appendix [A| we may now state and prove our main results 
on the existence, uniqueness and completeness of the geodesics in the regularised spacetime 

(pl), (M: 


Theorem 3.1 (Existence and uniqueness). Consider the geodesic equations (2.17) with initial 
dat a (|3.9[ ). Then for all e small enough (more precisely for all £ < £ 0 , where eo is constrained 
by (A.14)j, there exists a unique smooth solution = {Ue,Ve, Zpf) on [as,o:e + 1 ]], where p 
is independent of £ (and explicitly given by ( A.12 [) ). 



















GEODESICS IN NONEXPANDING IMPULSIVE GRAVITATIONAL WAVES WITH A, PART I 


9 



Figure 1. The ?7-component of the ‘seed geodesic’ 7 is depicted in black 
until it reaches the regularisation sandwich at parameter time t = i.e., 

U{as) = —£■ While in the background spacetime it would continue as the 
dotted red line to ?7 = 0 at t = 0, in the regularised spacetime it continues 


Theorem 
at t = /?( 


as 7 £ of (3.8) (depicted in green) solving the equations (2.17) with data (3.9). 

3.2 guarantees that 7 ^ (for e small) leaves the regularisation sandwich 
and continues as background geodesic 7 ^ of (3.14) with data (3.13). 


Proof. As noted in the appendix it suffices to first solve the (simplified) model system (A.l) 
for (ugjZg). Identifying {Us,Zp^) with (uejZe) the initial data (3.3) and (3.9) transfer to the 
data ( A.3[ ), (A.2). Then (3.10) implies ( A.4[ ) and Theorem A .6 applies to provide a unique 
smooth solution (I/^, Zp^) of (2.17), (3.9) on [«£,«£ + 7 ], with rj given by (A.12). 

Finally we solve the linear equation for I 4 to obtain the claimed smooth solution = 
{Ue, Ve, Zps) on [«£, ae + rj]. □ 

Next we make sure that the solutions just obtained, which by construction enter the wave zone 
at Ue = —£ at parameter time t = Ue with positive speed U^, in fact leave the sandwieh region, 
that is they reach Us = £ within their time of existence 7 . Consequently, they naturally extend 
to the background (anti-)de Sitter universe ‘behind’ the sandwich region. Observe that here 
it is vital that 7 in (A.12) is independent of e. 


Theorem 3.2 (Extension of geodesics). The unique smooth geodesics 7 ^ of Theorem 3.1 


extend to geodesics of the background (anti-)de Sitter spacetime ‘behind’ the sandwich wave 


zone. 
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Proof. Let 7 £ = {Us,Vs, Zp^) be the unique solution of (2.17), (3.9) given by Theorem 3.1 
By the definition of the ‘solution space’ (see ( A.5[)) we obtain 


(3.11) 


Us{0ie + 7]) = -£ + 


L 


a^+Tj 


U£{s) ds > —£ + > —£ + Se > £. 


since e < rjU^/Q by (A.13). 

So for such £ the geodesic ye leaves the wave zone and extends to a geodesic of the background 


spacetime since the geodesic equations (2.17) coincide with the geodesic equations of the 
background (anti-)de Sitter spacetime for Ue > £. □ 


Recall that by construction the global geodesic = {Us, Vs, Zps) with data (3.9) of Theorem 


3.2 for t < as, i.e. ‘in front’ of the sandwich, coincide for all (small) e with the single ‘seed 


geodesic’ y with data (3.3). However, ‘behind’ the sandwich the geodesics y^ for each e will 


coincide with a different geodesic of the background spacetime. To make this observation 
more precise, define the affine parameter time when the geodesic yg leaves the sandwich wave 
zone by jds. 


(3.12) Usit = /3s) = e. 
and denote the corresponding values of y^ at (3s by 

(3.13) yg(/3g) =: (e, ^0+, Z°+), %{(3s) =: (f7°+, 14°+, 

Then for t > (3s the geodesic y^ will coincide with the geodesic 

(3.14) y4 = (t/4,H+,Z+) 


7O+ 

^pe 


)• 


of the background (anti-)de Sitter space with the data (3.13), see Figure 1 and also Figure 2. 
Observe that the data (3.13) is normalised and constrained, more precisely we have: 


Remark 3.3 (Preservation of constraints and normalisation). The fact that the data (3.3) of 


the ‘seed geodesic’ y in (3.1) is constrained and normalised, i.e., it satisfies (3.4) and (3.5), 


implies that also the data (3.9) of y^ is constrained and normalised. Clearly these conditions 


are propagated by y^ being a solution to (2.17). Moreover, at t = /3g the regularised metric 
gs and the background metric g coincide and so the data (3.13) of y^ is constrained and 
normalised with respect to the background spacetime. 

While the preservation of the constraints confirms the consistency of our construction, the 
preservation of the normalisation, in particular, implies that the causal character of y^ (and 
y4) is the same as the one of the ‘seed’ y. 

Note that the geodesic y^ ‘behind’ the regularisation sandwich depends on e (only) via this 
initial data. Interestingly, as we will detail in Section for t > 0 the family of geodesics y^ 
of the regularised spacetime converges for e —)> 0 to a unique geodesic y+ in the background 


with data given by the limits of (3.13). This will explicitly describe the effect of the impulsive 


gravitational wave on the geodesics in (anti-)de Sitter universe. 

In the remainder of this section we will formulate completeness results for the regularised 
spacetimes. First we remark that actually our results allow us to make the smallness as¬ 
sumption on e precise: e has to be smaller than eo constrained by (A. 14). This, however. 


means that the specific e from which on a certain geodesic yg becomes complete depends on 


its data (3.3), i.e., on the ’seed geodesic’, and there is in general no global £ from which on 


all geodesics hence the spacetime is complete. A ‘global’ completeness result in the spirit of 
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Figure 2. The Z-components of two solutions 7 ^^ (purple) and 7^2 (green) 
of the regularised equations (2.17) with the same ‘seed geodesic’ 7 of (3.1) are 
depicted for ei > £ 2 - The regularisation sandwich is given by and 

respectively. The dotted red line represents the Z-components of 7 , 
while the black dotted lines are said components of 7 ^. 


however, can be obtained using the geometric theory of generalised functions ([HI Ch. 
3]) and we reserve its detailed presentation for future work. 

To formulate completeness results in our current setting, the dependence of e on the data 
discussed above also makes it necessary to be careful about global effects. Indeed, geodesics 
in the background spacetimes with ae > 0, that is spacelike geodesics in de Sitter space 
and timelike geodesics in anti-de Sitter space are periodic. Consequently geodesics 7 ^ in 
the regularised spacetimes constructed from such ‘seed geodesics’ 7 will share their causal 
character (see Remark 3.3) and hence cross the wave zone infinitely often and we (have to) 
(re)apply Theorems 3.1, 3.2 again and again. However, note that the geodesics may enter 
the regularisation sandwich region each time with different data. So the £ from which on 
Theorem 3.2 applies may in principal become smaller and smaller on successive crossings 
with no positive infimum. In such a case, given an initial geodesic 7 as in (3.1) and given 
any finite number N of crossings we can specify an £ from which on the geodesic 7 ^ extends 
to cross the wave zone A^-times. However, we cannot give a (global) e for which the geodesic 
7 e extends to all (positive) values of its affine parameter. Consequently we prefer to avoid 
multiple crossings of the impulse in the formulation of our results by restricting to causal 
geodesics in the de Sitter-case (neglecting unphysical tachyons only) and working with the 
universal covering spacetime in case of anti-de Sitter space. 

Note also that in our discussion so far (see the specification of the ‘seed geodesics’ 7 at 
the beginning of this section) we have exclusively dealt with geodesics with non-constant 17- 
component. Hence it remains to deal with geodesics travelling parallel to the surface {U = 0}. 
In case U = const 7 ^ 0 the geodesic will never enter the sandwich region of the regularised 
spacetime, provided e is small enough. Staying entirely on the constant curvature background 
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such a geodesic clearly is complete. To discuss the geodesics with U = 0, observe that the 
surface {?7 = 0} is totally geodesic (not only in the background but also) in the regularised 
spacetime, which can be seen from the ?7-component of the geodesic equations (2.17) (see also 
[32], the discussion prior to Thm. 3.6). Hence such geodesics have trivial [/-components and 
consequently the system (2.17) reduces to the background geodesic equations which again 
leads to completeness. So we hnally arrive at: 

Theorem 3.4 (Causal completeness for positive A). Every causal geodesic in the entire class 
of regularised nonexpanding impulsive gravitational waves propagating in de Sitter universe 
(i.e., (2.9), ( | 2 . 10 ) with A > 0 and a smooth profile function H) is complete, provided the 
regularisation parameter e is chosen small enough. 

Theorem 3.5 (Completeness for negative A). Every geodesic in the entire class of regularised 
nonexpanding impulsive gravitational waves propagating in the universal cover of anti-de Sitter 
universe (i.e., (2.9), ( 2 . 10 | ), with A. < 0 and a smooth profile function H) is complete, provided 
the regularisation parameter e is chosen small enough. 


Remark 3.6 (Non-smooth profiles H). In case the profile function H in the metric (2.9) is 
non-smooth — which, in fact, occurs in physically interesting models where H possesses poles 
on the wave front at z = ±1, see (|2.4|) — our method still applies but some care is needed. 


Indeed, if a ‘seed geodesic’ 7 hits the surface at [/ = 0 away from the poles we may work on 
an open subset of the spacetime with the poles of H removed. We only have to choose the 
constant Ci in (A.5) so small that the curves in the ‘solution space’ stay away from the 


poles, and then Theorems 3.1 and 3.2 still apply. The only ‘seed geodesics’ 7 which do not 


allow for such a treatment are those which directly head at the poles (i.e., 7 of (3.1) hits the 
pole at 7 ( 0 )), which is in complete agreement with physical expectations. 


4. Boundedness properties of geodesics 


In this section we establish boundedness properties of the global geodesics in the regularised 
spacetimes obtained in the previous section. In particular, we will prove local boundedness of 
7 e and of some components of its velocity uniformly in e. These properties will be essential 
in the next section where we derive the limits of 

To begin with observe that the fixed point argument of the appendix already gives uniform 
boundedness of the U- and Zp-components together with their first order derivatives. On the 
other hand, the H-component was not involved in the fixed point argument and we have to 
establish its boundedness properties using the H-component of the geodesic equation (2.17). 
Since this equation involves a A-term which is not multiplied by U^, the H-component of 7 ^ 
is not uniformly bounded in the regularisation sandwich. However, we will show that Vfifie), 
i.e., the H-speed when the geodesic leaves the regularisation sandwich is uniformly bounded. 


Proposition 4.1 (Uniform boundedness of geodesics). The global geodesics 7 ^ = ([ 4 , 14 , Zpe) 
of Theorem \3.i^ satisfy 

(i) [4 and Us are locally uniformly bounded in e, 

(ii) Zps and Zps are locally uniformly bounded in e, 

(iii) 14 is locally uniformly bounded in e, and 

(iv) Vsifie) is uniformly bounded in e. 

Observe that by Lemma |A.2| the time fis when the geodesic leaves the regularisation strip 
satisfies fis < «£ + 4e/[7° —)• 0. 
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Proof. Items (i) and (ii) are immediate from Theorem A.6 
To deal with the F-component we write 


t S 


(4.1) 


\V,{t)\ < |y0| + |y0| \t-ae\+ I I \V,{r)\ drds 


and estimate each term in the differential equation (2.17) for V^. To begin with we claim that 
(4.2) Vs is bounded on [as,/3s] ■ 


Indeed, proceeding as in the proof of Proposition A.3 and, in particular, using (A.5) and 
Lemmata A. [^we obtain 

t S 

\Vs{r)\ drds < Q||i7||^l||p'||^^(f70)2 + 3\\DH\\^-Jp\\^{\Z^\ + C 2 ) ^17°) 


2 /4e\2 


PKi/oJ 


9,.- 


+ ( 1 + piU^?\\Gs\\oc + '^UmH6sU, 


3.V 


\H\ 


00 llPIloo 


(4.3) 


9. ■ 


+ ^ ( 1 + ^(C/°)2||G,||oo + f-U^\\{H6sU, 


3 ■ 


t S 


8 


t S 


I Vs (r) I drds 


<C + C / \Vsir)\ 1 + 


Xeir) 


drds. 


where C is some generic constant and Xe is the characteristic function of [as, 13s]. Moreover, 
we have used that HGgHoo = 0(l/e) = [[{H6sUs)[[oo- So overall we obtain by a generalization 
of Gronwall’s inequality due to Bykov [3l Thm.ll.l]^ 

S 

(4.4) 


Vs[ < C{1 + e) exp ^ f f C[l + -)drds^ < C < 




establishing the claim ( |4.2[ ). 

Now we may prove (iv): Writing Vs{/3s) also as an integral and using boundedness of Vs on 
[as,f3s] we may proceed as in ( |4.3| ). Then again we obtain uniform boundedness of all the 
terms but the first one, which involves the d'-term. To estimate this one we use integration 
by parts to obtain 

j\{Zs{s))5',{Us{s))Ul{s)ds = j H{Zs{U-\r)))6',{r)UsiU-\r))dr 

— £ 

£ 

(4.5) = 6s{r)H{Zs{U-Hr))) UsiU-\r))\[^ - J 6s{r)(^H{Zs{U-\r))) Us{U-\r))) dr. 


^Note that it suffices that Xe is integrable rather than continuous. 
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Now the first term vanishes and the second one is bounded independently of e by (i), (ii) 


the fact that ^ is uniformly bounded away from zero by (A.5), and since C/g is uniformly 


bounded by (2.17). This establishes (iv) 


To prove (iii) it remains to show that 14 is bounded on any compact interval disjoint from 
(«£, j3e)- But this follows from the fact that ye outside of (oe, I3e) solves the geodesic equation 
of the background spacetime of constant curvature, and that 7e(/3e) and jeiPe) are uniformly 
bounded by (i), (ii), (iv) and (4.2). □ 


5. Limiting geodesics 

In this final section we consider the limit e —)• 0 of the unique global smooth geodesics % 
of the regularised spacetime (2.9), (2.10) obtained in Section(Thm. 3.1 


Thm. 3.2). This 


physically amounts to explicitly determine the geodesics of the distributional form ( 2 . 1 ), ( 2 . 2 ) 


of all nonexpanding impulse gravitational waves propagating in (anti-)de Sitter universe. 
In particular, we will prove that the geodesics 7 e converge to geodesics of the background 
(anti-)de Sitter spacetime with appropriate but different data on either side of the impulse 
(U < 0 and U > 0 respectively), which from a global point of view amounts to convergence of 
7 e to geodesics of the background which have to be matched appropriately across the impulse. 
The technical calculation of the limits is given in Appendix [B| 

To make our claims on the convergence of 7 e = {Ue,Ve, Zp^) precise we introduce the following 
notation for the geodesics of the background (anti-)de Sitter universe: Let 7 = (17, V, Zp) be 


a ‘seed geodesic’ as in (3.1), that is 17(t = 0) = 0 and 7 assumes the data (3.3), i.e., 

(5.1) 7 ( 0 ) = (0,yO,ZO), and 7 ( 0 ) = (17° > 0,1/°, 4°), 

41) and are normalised as in (3.5). Furthermore 


7(0) = (0,y^4), 

where the constants satisfy the constraints 

let y"'' = {U~^,V~^, Z^) be a geodesic of the background again crossing U = 0 at t = 0, i.e., 
with U~^{t = 0 ) = 0 and with data 

(5.2) 7+(0) = (0,B,zO), 7+(0) = (170,C,Ap) 

where we define 


(5.3) 


Ap := limZpeifde 
£-:>0 


B = lim Vs{j3e ), 
£-^•0 


C = lim V^{j3e ). 
£-^•0 


Recall that (3^ < + As/U^ —>■ 0 is defined to be the time when the regularised geodesic ye 

leaves the regularisation strip, i.e., 17e(/3e) = s. Finally define 7 = (17, V, Zp) by 


(5.4) 


7 ( 1 ) : = 



t < 0 
1 > 0 . 


We will show that ye converge to the ‘matched geodesics’ 7 of the impulsive spacetime, which 
from now on we will also call ‘limiting geodesics’ with ‘past branch’ 7 and ‘future branch’ 
y"*", see also Figure 3. Note that the respective notion of convergence of the individual 
components of ye will differ, subject to the regularity of the respective components of the 
‘limiting geodesics’. Indeed, 7 = {U,V,Zp) has a smooth first component 17, while Zp is 
continuous with a finite jump in Zp (determined by Ap) across the impulse at 1 = 0, and V is 
even discontinuous across 1 = 0 with a finite jump in V and V (determined by the coefficients 
B and C, respectively). 


Observe that at the moment we only know the limits in (5.3) to exist for subsequences (by 


uniform boundedness of Zp^, I 4 and I4(/3e), cf. Proposition 4.1) and hence Ap, B, and C need 
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not be uniquely defined. We will, however, prove convergence and we will derive an explicit 
expressions for Ap, B and C in Proposition 5.3 below. But first we state and prove the main 
assertion on the limits of the geodesics in the regularised spacetime: 

Theorem 5.1. The geodesics ye = {U£,Vs, Zp^) of the regularised spacetime derived in The¬ 
orem's^ converge to the ‘limiting geodesics’ ^ of (5.4) in the following sense: 

(i) Ue^ty ^nC^ 

(ii) Zp£ —>■ Zp locally uniformly, 

Zp in distributions and uniformly on compact intervals not containing t = 0, 




pe 


(iii) —)■ V in distributions and in on compact intervals not containing t = 0. 

Observe that Zp is discontinuous across t = 0, hence convergence of Zp^ cannot be uniform 
on any interval containing t = clh and the same holds true for and Vs. 


Proof. First we consider the U-component of on the interval [ae,/?^], where we have from 
the geodesic equations (2.17) resp. (2.6) 


(5.5) 


\Us{t)-U{t)\<e+\e\ 


+ 



Us 


u 


p - UfHds 

rV 2 


dr ds 


Wsiis Gs - Usiis inSsUsY 


a a? - UfHSs 


To estimate I observe that (cf. the proof of Lemma A.l) 

1 


(5.6) 


aa 


2 - mH5s 


1 


cra^ 


< -,e\m 


oo IIPI OO 


dr ds =: e + lell + II. 


< Ce, 


with C a generic constant, and consequently by (A.5) 

rt rs 


I < 


Us 


Us 


a a? - UfHds a a? 


(5.7) 


< rYis + Ci)Cs H —K 


S 

t rs 


+ 


Us 


u 


Us-U 


aaf (To2 


dr ds. 


dr ds 


For the term II we obtain, as in the proof of Proposition A.3 (cf. (A.16), (A.18)), that 
II < Ce. Hence, overall 


(5.8) 


\Us{t)-U{t)\ <C£+^ 


Us-U 


dr ds, 


and so again by Bykov’s inequality \Us{t) — U{t)\ = 0{s). In the same way we see that also 
\Us{t) — U{t)\ = 0(e) and so 


(5.9) 


sup \Us(t) - U(t)\ + \Us(t) - U(t)\ ^ 0 . 


We now turn to the Zp-components on the interval [ae,/?^]. We have 
(5.10) 

sup \Zps(t) - Zp(t)\ < sup \Zps(t) - Z° I + sup |Z° - Zp(t)\ -^0 (e 0), 

as<t<fiE as<t<fiE as<t<j3s 


Wnless t = 0 is the right endpoint of the interval. This, however, is rather an artefact due to our choice 
setting 7(0) = 7(0), cf. (5.4l. 
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where = Zp{ae) was defined in (3.9). Indeed by continuity of Zp the second term converges 
to zero since Og, /?£ —)• 0 and the first term can be estimated using the differential equation 
by 


\ZpS)-zU < 


rl^e fS 

J J 

J OLfr J OLp 


DH5eU, 


eZl + \ZI,U, Ge - ZIU, {HdeUeY 


aa 


2 - U^H6, 


dr ds 


(5.11) 


+ \ Zpe\{c>-e + 13s) < Ce . 


Here we have used that the inner integral is bounded by (A.27)-(A.29). 

Now we finish the proof of (i) and establish the claim on uniform convergence in (ii) and 
(iii). First note that by construction there is nothing to show for t < 0. For t > /Sg we use 
continuous dependence of solutions to ODEs on the data. Indeed for such t both 7 ^ and 
7 = 7 “*' are solutions to the same differential equation, however, with different data which is 
given for 7 ^ at t = /?£ by (5.2) and for 7 + at f = 0 by (3.13). More precisely, for all T > 0 
(which implies T > (3^ for e small) we have 


sup {\'yg{t)-'y+{t)\,\%{t) - 7 +(t)|) 

Pe<t<T 


(5.12) 


< max [\-fe{l5e) - 7“^ (/?£)!, \is{l3e) - 7^(/3e)l) e' 


TL 


where L is a Lipschitz constant for the right hand side of the geodesic equation of the back¬ 
ground on a suitable compact set. Note that such a set exists by the boundedness properties 
of 7 £ established in Proposition 4.1, i.e., YeWs) is uniformly bounded. Finally, for the terms 
in the maximum in ( |5.12 ) we have 

(5.13) \Ug{Yie)-U+{^g)\^Q, \Ug{^g)-U+{Yg)\^^ by(l5;§, and 

\Zpe{l3e) - (/?£)! 0 by ( |5.10[ ) , 

whereas for the remaining terms we write 


(5.14) 


\Vg{l5g) - V+{l5g)\ <\Vg{Pg) - H+(0)| + |H+(0) - V+{l5g)\ , 
\Vg{l3g) - 1/+(/3£)| < \Vgm - H+(0)| + |H+(0) - V+{l3g)\ , 
\Zpg{^g) - zt{i3g)\ < \Zpg{^e) - ztm + \zt{Q) - 4+(/3.)| . 


Now in each line the last term on the right hand side goes to zero by smoothness of 7 '*', while 
for the respective first terms we have by our choice of data (5.2), (5.3) 

|F£(/3£)-1/+(0)| = |K(/3£)-B|^0, 

(5.15) \Vg{l3g) - 1/+(0)| = \Vg{Yis) - C| ^ 0, 

\Zps{l3s) — Zp (0)1 = \Zpg{tg) — Apl —)■ 0 . 


Finally to prove the distributional convergence in (iii) thanks to the uniform convergence of 
Vg established above we only have to consider the integral 
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for a test function (p on M. This, however, converges to zero by the local uniform boundedness 
of 14 established in Proposition |4.1[ iii). In case of the distributional convergence in (ii) we 
argue precisely in the same manner, now using Proposition |4.1[ii). □ 


Remark 5.2 (Normalisation and constraints in the limit). The convergence result provided by 
Theorem |5 .1 1 also guarantees the preservation of the normalisation, i.e., the normalisation of 


the ‘seed geodesic’ 7 , which by Remark 3.3 carries over to the regularised geodesics 75 (and 
hence to 7 ^), also carries over to the ‘future branch’ of the ‘limiting geodesic’ 7 ^. To see this 
we just write 

e = ge{le{Pe)){ie{Pe),ieWe)) = g{ieWe),ie{(ie)) 

(5.17) = -2UM)Vem + {Z2e{(ie)f + + ^{ZieiHe)? 

^ -2I7°C + A22 + A32 + cjA42 = 5(7+(0),7+(0)). 


Here the first equality follows from Remark 3.3 and the second one follows from the fact, 
that at the regularised metric agrees with the (constant) background metric. Finally, 

convergence is due to Thm. |5.1[ i) and our choice of the data ( |5.2[ ), (5.3). 

Also by a similar (actually simpler) argument the constraints carry over from 7 ^ to , which 
again confirms consistency of our construction. 


To end this section and the entire paper we now explicitly evaluate the limits in (5.3), thereby 


showing that the ‘limiting geodesics’ (5.4) of the smooth global geodesics 7 ^ of Theorem 3.2 


the regularised spacetime (2.9), (2.10) coincide with the geodesics (2.6) of the distributional 


spacetime ( 2 . 1 ), ( 2 . 2 ) derived previously in m- 


m 



Figure 3. The P-components of 7 ^^ (purple) and 7^2 (green) for si > £2 are 
depicted. They converge to the the ‘limiting geodesic’ 7 whose ‘future branch’ 
is separated from its ‘past branch’ 7 (black outside and dotted red inside 


the regularisation sandwich) by the ‘jump’ B calculated in Proposition 5.3 


The ‘jump’ C in P is indicated by the different P-slopes of the ‘past branch’ 
7 and the ‘future branch’ . 
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Proposition 5.3. The constants Ap, B, and C determining the data for the ‘future branch’ 
7 + of the ‘limiting geodesics ’ 7 are explicitly given by 




(5.18) 


B = YmVM) = Ih{Z^^) + V\ 


£—^0 




0^2 


C := lim VM) = y +^{ H,2{Z^rY + + aH^Zf!) 

£—^-O 5 \ 


+ ^F(Z0)2 - (6P‘^Z^IIjZ^)f) 

aa^ aa^ ^ ^ ' j 


t/' 


^ - i7(ZO)) po + ^dPm^.iZ'^YZ^ , 


The calculation is rather technical and we sketch the main points in Appendix]^ Finally 
remark that our results are fully compatible with the ones in m- 


we 


Remark 5.4. To simplify the comparison of the results on the ‘limiting geodesics’ of Propo¬ 


sition 5.3 with the heuristically derived geodesics of the impulsive wave spacetime of m 
eqs. (38),(39)] we remark that there the geodesics were restricted to = Z^ = 0, and 
= 1. Moreover recalling that l/{aa^) = A/3 and using the notations G(0) = G(Zp) = 


5P^Z^H^g{Z^) - H{Z^) and Fr( 0 ) = H{Z^) of [HI, equations (|5.18[) take the form 


Ai= Vf,(0 )-|z/G(0 ) 


k, = UaH^M-^ZlG{t)) 


B= 2^(0)> 


(5.19) C= ^(^,2(0)2+//,3(0)2 + uF4(0)2 + |i/(0)2-^(5^’%Oi/,,(0))') ^ 

Finally taking into account that e = sign(A) in [21] as well as the slightly different definition 


of C in |3ll eqs. (37)] we see that (5.19) indeed agrees with eqs. (38) of [31] 


Summary 

In this paper we have rigorously investigated all geodesics in the entire class of nonexpanding 
impulsive gravitational waves propagating in (anti-)de Sitter universe, extending thus many 
previous studies of geodesic motion in the class of impulsive pp-wave spacetimes with vanishing 
cosmological constant. Following m we employed the distributional form of the metric in the 
context of a b-dimensional embedding formalism. We have applied a regularisation technique, 
replacing the Dirac-(5 by a general class of smooth functions—the model delta nets (2.8). Since 


we have never used any special property of the regularising net, our results are completely 
regularisation independent within this class. In physical terms this means that the formal 


distributional form of the impulsive metric ( 2 . 1 ), ( 2 . 2 ) is understood as a limit of a family of 


spacetimes with ever shorter but ever stronger sandwich gravitational waves (2.9), (2.10) of 
an arbitrary smooth profile 6^. 
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Although the resulting regularized geodesic equations (2.17) form a highly complicated cou¬ 
pled system, we were able to prove in Section the existence and uniqueness of geodesics 
crossing the wave impulse, leading to completeness results, see Theorems 3.1, 3.2 and 3.4, 3.5[ 
Observe that, in particular, we prove that the geodesics of the regularised spacetime hitting 
the wave zone and hence interacting nonlinearly with the impulse actually cross it. This is a 
physical information not provided by the approach of m in which this feature is built into 
the heuristic ansatz. 

Our proof is based on the application of a fixed point theorem, the technical details of which 
can be found in Appendix]^ There we have extended the range of applicability of this kind 
of fixed point techniques, originating in |22j and generalised in [36], to a far more involved 
situation where the higher order ‘contraction estimate’ contains a term proportional to 1/e, 
cf. (A.32). In this way we have pushed on the crucial point of the method to the estimate 
(A.29), cf. the remark in the middle of page |25| below the proof of Proposition |A.3| This raises 
hopes that these techniques can be extended to the even wilder ‘4-dimensional’ distributional 
form of the metric (1.6) in the future. 

In Section we studied boundedness properties of the global geodesics in the regularized 
spacetimes. These technical results, summarized in Proposition 4.1, were essential for per¬ 
forming their limits in the final Section Due to the complexity of the system of geodesic 
equations (2.17) it seemed advisable to simplify the ‘usual arguments’ in this limiting pro¬ 
cedure. We have done so by abstracting from the concrete form of the (limiting) geodesics 
and repeatedly using continuous dependence of solutions of ODEs on its data. In this way 
we were able to show that, as the regularisation parameter goes to zero, the solutions of 
the regularised geodesic equation converge to unique geodesics of the background (anti-)de 
Sitter spacetime (Theorem 5.1) which have to be matched appropriately across the impulse. 
In fact, we rigorously derived the explicit form of these matching conditions (some overly 
technical related calculations are contained in Appendix]^. The resulting coefficients ( |5.18| ) 
of Proposition 5.3 fully agree with previous results derived by a heuristic approach in |31j. 
Remarkably, the impulsive limit is completely independent of the specific regularization, i.e., 
in the limit e —)• 0 it is the same for any smooth profile of the sandwich gravitational waves. 
Finally, as mentioned in the introduction in [32], we have recently investigated the complete 
family of nonexpanding impulsive gravitational waves propagating in spaces of constant curva¬ 
ture (Minkowski, de Sitter and anti-de Sitter universes) employing the (Lipschitz) continuous 
form of the metric (1.4). Using Filippov’s solution concept for differential equations with 
discontinuous right hand side we proved existence and uniqueness of continuously differen¬ 
tiable geodesics. In section 4 of [32] we explicitly derived such geodesics using a C^-matching 
procedure resulting in specific matching conditions, namely equations (4.4)-(4.10) of [32j . 

A natural question thus arises about the mutual consistency of the two results, both obtained 
in a rigorous way but starting from two different forms of the metric, namely the continuous 
form of the metric (employed in [32]) and the distributional form of the metric in the context 
of the 5-dimensional embedding formalism (employed here and in |31ji. In fact, it was shown 
in the recent work [20] that the matching conditions of |32j and |31j are fully equivalent 
when appropriate coordinate transformations are applied. This result confirms that both our 
approaches are consistent. It follows that the understanding of geodesics in the complete 
family of spacetimes with nonexpanding impulsive gravitational waves and any cosmological 
constant now rests on firm mathematical grounds. 

These results set the stage for a sound mathematical analysis of the ‘discontinuous coordinate 
transformations’ between the continuous and the distributional forms of the metric. Together 
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with the results of [20] is seems now feasible to rigorously relate the continuous form of the 
metric (1.4) to the ‘5-dimensional’ distributional form (2.1), (2.2). On the other hand the 
technical advances of the fixed point techniques made here might eventually bring into reach 
a direct approach on the mathematical intricacies of the transformation (1.5). 
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Appendix A. The fixed point argument 


In this appendix we detail the fixed point argument used to prove a suitable existence and 
uniqueness result for solutions of the regularised geodesic equations (2.17) with data (3.9) that 
additionally guarantees the solutions to live long enough to leave the regularisation sandwich. 
To do so, we only have to prove existence of the C4 and Zp^-components, since the equation 
for 14 decouples and is linear, hence can then be solved on the domain of existence of {Us, 
Zpe)- Moreover, the sign-difference between the Zj^-equations and the Z 4 £-equation can safely 
be ignored in the estimates leading to the fixed point argument. Therefore, in this appendix, 
we only (have to) deal with the following simplified model system 


(A.l) 


Ue 


A 


1 

2 


DHdsiil 


-{e + ]j^ulGs-Ue {H Ss 
-{e + ^ujGs - tie {H 6s 


Ue 

(Ta2 - ulH6s ’ 

^e 

(Ta 2 - ulHds ’ 


where H = H{zs) is a smooth function on DH denotes its gradient, and Gs{us,Zs) := 
DH{zs) 6s{us) Zs -|- H{zs) 6^(us) Us- We will also frequently use the notation Xs = {us, Zs). 


We begin by setting up the initial data. Let rj > 0 and let Js = [o-s, o-s -U rj\ be the parameter 
interval where we look for solutions. In accordance with the strategy employed in Section]^ 
we will pose initial data at t and compare it to fixed data (corresponding to the initial 
data of the ‘seed geodesic’ at t = 0). So let 

= (u°, G M X and 7:° := (t6°, i°) G M x be given and set 
(A.2) Xs{as) = {us{as),Zs{as)) = {ul,z^^) and Xs{as) = {ue{as), Zs{as)) = {ul, z^) 


and let additionally 

(A. 3) u^,u^ G M, z^, 7° G be given and write := {vP, z^), x^ := (tt°, 7°) . 

As detailed in Section we exclusively deal with data satisfying 


Ug = —e, ttg > 0 and = 0, ir > 0, with the additional assumption (cf. (3.10)) 
(A.4) x^ and as e —^ 0. 
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We will apply our fixed point argument on a complete metric space which we will call the 
‘solution space’ and which is given as the closed subset of defined by 


Xe ■■= jxe = {Ue,Ze) G C^(j£,]R^) : Xeia^) = Xeiae) = and 


(A.5) 


Xp — X 


— ^1) \\^£ ^ ||oo ^ ^2; ^ 


rl 


.-.o 


-U , -u 


.2 ’ 2 


}■ 


Observe that we have ‘centred’ the functions in Xp around the ‘fixed’ initial data (A.3), while 
the prospective solutions are required to assume the e-dependent data (A. 2 ) at t = ap. Also 
note that the final condition forces iip to stay positive, which is the essential ingredient that 
forces the solutions to leave the regularisation sandwich. Now we arrange the constants as 
follows: First let Ci > 0 and set 

C2 :=l + maxLu^\\DH\\oo\\p\\i, + <^i) (311^^^110011/91100(1^°! + C-i) + ll^lloo||p'||oo) , 


(A. 6 ) 


^(|z°| +Ci) ( l + ^u°||^ 


P 


+ 


where ||^/^||oo and \\DH\\oo are taken over the closed Euclidean ball Bc^{z^). Also p is as in 
(2.8). Observe that the space Xp only depends on e via the domain and the initial data 

jg. 

Next we define the solution operator Ap acting on Xp for all t € via Ap{xe){t) := {Al(xe){t), 
A'^{xe){t)) with 

(A.7) 


Al{xe){t) := - 


A^^{xe){t) := 


ff 

J Ois: 'J (Xs: 

ct ps / 


eUs + \ueulGe - UpUe {H5eUe)' 


aa^ — uiH 5e 


dr ds + u^^{t — ap) — £ , 



T~\ TT r • 2 ^ 2 

-DHOeU^ -^ 


eZe -I- \zeulGe - Z^ile {HdpUe)' 


a a? — u^HSe 


dr ds + i°(t — Up) + z° , 


where again we have suppressed the dependence of 5^, G^ and H as well as their derivatives 
on the variables. 

Our first step will be to show that the operator Ap takes X^ to Xp , see Proposition |A.3| below. 
We begin with two preliminary results. First we bound the term in the denominator of Ap 
from below. 


Lemma A.l. Suppose z G 5c'i(^°) then for all u G 


(A. 8 ) 

provided e < a^/( 2 ||/ 9 ||oo||F 7 | 


1 2 
aa? — H{z) 5s{u) ~ ’ 


Proof. First, in case jttj > e we have u 4 supp(5e) and consequently -j— 

\aa^ — u^HOp\ 

Second, in case juj < e we have \u^H{z)5s{u)\ < e||i7||oo||/0||oo < a^/2 and therefore in both 
cases 

1 2 

jcra^ — u^H5e\ ~ ' 

□ 
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The second preliminary result shows that the conditions on imposed in (A.5), i.e., that 
Us > h^/2, prevents the u-component from slowing down too much in the sense that Us{t) 
leaves the sandwich region early enough. To state the result in a precise way we dehne for 
Xe = (rig, Zg) G Xg the set 

(A.9) rg(xg) = rg(Rg) ;= {t e Je |rr£(t)| <e}CJe, 


which is the maximal set where the terms in (A.l), (A.7) involving 5^ or 6'^ are non vanishing. 
We now have: 

Lemma A.2. The diameter ofT^{uc) is bounded for all Xg = {us,Ze) G 3£g by 


(A.IO) 


diam(rg(rig)) < - 


4e 


Proof. For Xg G Xg let t G rg(xg) which implies \u^{t)\ < £ and so 


(ATI) 


1 


e > rtg(t) = tig + / Ug(T) dr > Mg + -ir{t - Og) 


But this implies t < Og + 2 (e — u^)/vP = ccg + Ae/vP. □ 

Now we may state and prove that Tg(Xg) C Xg provided rj is chosen appropriately and e is 
small enough. 


Proposition A.3. Set 

_ . [ a2 Cl 2Ci a^Ci a^ C 2 

•“ “'""I ’ 24^ ’ |T^ ’ M\\p\\i\\DH\\ocU^ ’ 12 (|z 0 | + Cl) ’ 8 (|z 0 | + Ci) ’ 

^(m°(|.°| + Ci){3\\DH\\Moo{\z^\ + gi) + l|g||oo||p^||oo))~\ 

(A.12) 

Y2 + ^ 1 ) (3||T>77||oo||p||cx)(|i^°| + C 2 ) + -m°||77||oo(||/5'||cx) + IIt’IIcxd)^^ | 

and 

(3|l"»ll~ll'>ll~(l‘-"l + C.) + ll^flloollp'lloo) ■ , 

(A.13) ^ (^3||Z)F|U||p||oo(|iVC2) + 2 «°ll^lloo(||p'||oo + ||p||oo)) , 

{3\\DH\U\p\U\A+C2))~\'^,hY 


Now choose eo such that 
0 < £0 < ^0 ) 

(A. 14) |Mg — m'^I < , \z^ — z^\ < ^ , and \z^ — z^\ < 1 for all 0 < e < eo • 

8 6 

Then for all £ < £0 the operator Ag maps Xg to Xg. 


Observe that by (A.4) there exists eo, that guarantees the estimates in (A. 14) to hold. 
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Proof. We begin by estimating 

^ eUe + \ueulGe-UeUe{H5eUe)' .g 

(A.15) = - f - -do + “? 


and proceed term by term beginning with the latter two under the integral which we will see 
1 as e —)• 0. Indeed we ha^ 

Uei4{DH5sZe + H5'^Us) 


to vanish as e —)• 0. Indeed we have for all e < eo by the definition of (A.5) and by Lemma 

ED 


2 {aa? - ulH5e) 


-ds 


(A.16) 


< ^diam(r£(rte)) e ^SllDifUooll/olloo + Ci) + ||-fl||oo^||/o'||oo 


QoyO 

<^e{3\\DH\ 


cxD IIPMcxd \ 


3(|2°|+C'i) + ||i/||oo||/5^||oo) < g, 
where for the second inequality we have used Lemma A.2| and for the third that 


(A.17) 


eo < 


72vf 

Similarly, we have 

Ueils {H5eUeY 


{nDH\ 


OO P OO 


+ Cl) + ||i7||oo||/5^||oo) 


-1 


aa^ — u^H6e 


12 


ds 


(A.18) <^e(3£\\DH\\ 


loo IIP ICXD 


{\z^\+C2) + -u^\\H\ 


+ 


< o > 


where the final estimate again follows from our assumptions on £q. 

Finally to estimate the first term under the integral in (A.15) we write for Ue G Xg 

t 

(A.19) 


/ 3 3 

Ue{s) ds < —£ + T] -vP < -iPr] . 


Since —e < and by the last condition on eo in (A.14) we obtain |u£(t)| < and 


hence 

(A.20) 


eup 


2 - 


ds 


aa 


2 /■* , . M , 3 n 2 3 n 1 

< ^ / U£(s) ds < -^u V < - 2 ^ V < - 


where we have used that rj < 1 and ry < a‘^/{24u^), cf. (A. 12 ). Thus, by In^ — < | we 


obtain overall ||^Ag(x£) — u'noo 


.•, 0 | 


< i.e., ^Al{xs)it) G [hiP, |u°] for all t G J^. 


Moreover, using the above estimates, integrating once more and using e < ry we find that 
(A.21) ||Ae(a;£) - u°||oo < e + ry^ + ryw° < v{^ + u°) < Ci, 

due to the assumption that ry < Ci/(| + 
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Now we turn to the ‘spatial component’ of the solution operator. We have to show that 
\\Al{x,)-z«« 


(A.22) 


|oo 

rt rs 



-DHSeul - 

r) ^ c 


eZe + IzeiilGe - ZsU£{H5eUe)' 


a a? — u^H6e 


^ dr ds + z^{t — Ue) 


<Ci 


and again proceed term by term. To begin with we note the following auxiliary estimate 

rt r, rt A,0 


r 2 r ir 

(A.23) / Mu£{s))\ds =-J \6,{u,{s))\-ds 

Ja.fr ^ J a^ ^ 




\6^{u^{s))\usis)ds = 




rue{t) 2 

I < 




Now we have once more using the definition of 


(A.24) i 


ff 


DHS^Ug drds 


<-jjWlil|oJd|| 


V = ^ IIpIIi \\DH\\^u^t] < ^, 


where we have made use of r/ < Ci/((2/54)||/9||i||L)//||ooU°). Similarly since ?7 < {a^ Ci)/{12{\z^\ + 
Cl)) we obtain 


(A.25) 


t ps 



ezp 


ia,Ja, (JO? -ulH5e 


dr ds 


<|(|2»I + c,).,^<4(|/| + c.),< A. 


Furthermore, we estimate 


'a, J a, 0 -“^ - uiH 6e 


< 4h°(|z°|+<^i) mDH\ 


oo llPIloo 


{\z°\+C,)+ \\H\U\p Woo) ri<^, 


where we have used ry < +Ci)(3||OFf||oo||/o||oo(|2;°|+Ci) + ||Ff||oo||/o'||oo)) , and 

finally 


t ps 


Zsile {H5eUe)' 
(TO? -u‘lH5e 



dr ds 


(A.26) 


— +1^1) ^3e||Z)Fr||oo||p||oo(|^^| +C2) + -u^||Ff||oo^||p^||oo + IIpIIoo^^ P 


- 6 ’ 


where we have used the hnal condition on rj in ( A. 12 [ ). This establishes (A. 22 ) using the one 
before last condition on ry in (A. 12 ) together with \z^ — z^\ < Cij6. 

It remains to show ||^Ag(x£) — i°||oo < C 2 . As in ( |W24l ), (|A^ we estimate 
1 


(A.27) 


DH6^{u^)u^ ds 


< jIWi 


IDidlUti" < Y • 




aa^ — u^H 5 e 


ds 


<|(| 2 »|+C 0 d<^ 


where we have used the first condition on C 2 in (A. 6 ) and the sixth one on rj in (A. 12 ). 
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For the remaining two terms we have 

ds 


1 z^iilGe 


a a? -ulH5e 


(A.28) 


36 


[nDH\ 


oo llPIloo 


(|z°| + Cl) + ||Ff||oo||/5^||oo) — 


where we have used the second condition on C 2 in (A. 6 ), and 
Zeile {H6 sU£ 


f 

J OLp 


aa‘ 


— uiH5f 


ds 


(A.29) 


12 

- a 2 


{\z^\+C,)(3e\\DH\\ 


00 P 00 


(|i°l+C' 2 ) + 2 ^°||^|| 


+ 




00 M < 4 • 


Here we have used the fourth condition on eo in (A. 14) as well as the final condition on C 2 
in (|A. 6 |). □ 


Observe that in the estimate (A.29) it is absolutely vital that the term in the second line 
involving C 2 is proportional to e —otherwise we would end up in a circle and our method 
would fail. 

Our next step is to prove that the solution operator has a fixed point on Xg. To this end 
we need the following technical preparation. 

Lemma A.4. There exist constants C and C' (independent of e) such that for all X£,x* G 
we have 

rt 


(i) 


{5e{Ue)Ue - 4«X) ds 


< C||ue — u*||oo , and 


(ii) 


{Ki^sWe - KK)iKf) ds 


< C'llue - n*||oo ■ 


Proof. To prove (i) we first consider the case ||m£ — u^Hoo < £• We have by 


{SeiUe)u£ - 6£{u*)u*)ds 
(A.30) 


< 


f \6£{u£)u£ - d£{u£)ul\ ds + f \6£{u£)ul - SeiuDull ds 
J CXr J OLp 


2 , 
- hol 


L 


|(5e(Ue) - (5eX)| l^l ds . 


1111 y^£ 11 00 “t“ 

/r£(a:£)ur£(a;*) 

Now the last integral is non-vanishing only if 17/^1 < £ or |?/*| < £ hence we have in any case 
17/*I <26. Since both x* E Lemma A.2 applies so that using a mean value argument 
we may bound the integral by ^||/9'||oo||i^£ — ^^elloo- 


In case ||u£ — u*||oo > e we obtain again from Lemma A.2 


{6s{us)us - 6s{u*)ul)ds 


< 


[ |5£(n£)u£| ds + f \6 e{ul)ul\ds 

Jrpx,) 


irpx*p 

8 8 

(A.31) < ^ IIpIIoo £ < IIpIIoo llll-e ~ llcxo 

So we may chose C = ^max(||/ 3 ||i + 4||/9'||oo,4||/3||oo) and (i) is proved. 
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(h) is proved analogously with the choice C = max (4||/?'||oo + ||/9"||oo , 2||/9'||oo). 


□ 


We finally prove the key estimates which will allow the application of Weissinger’s fixed point 
theorem. 


Proposition A. 5. There exists a sequence of positive real numbers {an)n (depending on p, 
p', p", H, DH, D^H, and iP but independent of s) with XlnsN°° such that for all 
X£,x* G Xs with s < eo of (A.14) and p as in (A.12) and all n G N we have 


(A.32) ||(^.r(x.) - {A^nxDiy < - x:||ci . 

Proof It suffices to show ||Ae(xe) — A£(x*)||ci < {C/£)\\x£ — x^Wci for some appropriate 
constant C, since for higher powers we then may use 


(A.33) 



1 dtdti . . . dt2n-l < 


p 


2n 


(2n)! 


to obtain a converging series. 

We again proceed term by term, skipping some of the details of the (by now) routine estimates 
and only stress the technical key points. 

We start with the first term in WA^ipXs) — A^{x*^)\\qi. By writing the two summands on a 
common denominator we obtain 




__ds 

{u%YH(zl)5e{u%) 

\u£{ulfH{zl)5£{ul) - 
ululH{z£)5£{u£)\ds 


Ue{y*£fH{Ze)de{u*£)\ds 


(A.34) 

< \p\\u£ - U*||oo + \{\u^\ + Cif (lAp{H)\\z£ - <||oo^||p||oo + \\H\\ocC\\u£ - M*||oo) 

- + C'l)^(Lip(^) Jjll/5||oo + ||i?||oo(5)^ Ike - X*||ci , 


where Lip(i7) denotes the Lipschitz constant of H on Bc^{z^) and C is the constant given 
by Lemma |A.4[ 

For the second term we need the following auxiliary estimate which is proven by a combination 


of (i) and (ii) in Lemma A.4 


(A.35) 



G*ut\ds<C'\\x£-x;\\ci, 


where C' = \\D^H\U\p\\ooi\z^\ + Ci) + \\DH\\^ ((k0| + Ci)C + ||p||oo + ||p'||oo) + II^HooC'- 
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Abbreviating Cq := 3\\DH\\^ 

f \ueulG 


(Iz^l + Cl) + ll-H^llooll/s^lloo we are able to estimate 
ds 


\ul{ulfG* 


2 

< — 




aa? - u^H{ze)5e{ue) aa^ - {u*yH{z*)Se{u*) 

r*i 

^ UeuiG-u*{ul)'^G* + 


18(ttO)/a2G' a^Gg G' 


^ a* 


+ +til'll'”' 


X \\Xe - X^Wcl , 


by using (A.35) and Lemma A.4 


GG, 


G I 


H\ 


+ ^Lip(iL)||p||oo + -^\\H\\ 


ds 


The final term in i.e., 


(A.36) 


UeU. 


'^e(^H{Zs)6e{ 


Ue)Ue 


KK 


{Hiz:)5e{i 


aa? - u^H{ze)6e{ue) aa^ - {u*yH{z*)6eiu*) 


ds 


can be estimated in perfect analogy to the previous terms inserting and subtracting appro¬ 
priate terms wherever necessary to arrive at an estimate proportional to ||x£ — x*||ci. 

The ‘spatial component’ of the solution operator can be treated in a similar way. The 
only new aspect when estimating \\A^{xe) — A^{x*A\c^ is the following. When bounding terms 
like |G — G*| by multiples of ||x£ — x*||ci we hnd that they are no longer multiplied by 
and u*, respectively. Thus we cannot use the auxiliary result (A.35) and consequently terms 
proportional to 1/e remain. (Note, however, that the occurrence of 1/e-terms at this stage 
causes no problem at all in the application of the fixed point theorem, see below.) Summing 
up we arrive at 




< ^G\\xe-xl\\ci 


where G is some constant (as above depending on H, p, etc.). Furthermore, since 77 < 1 we 
obtain the same estimate for the zeroth order, i.e., \\AAxe) — Ae(a:*)||oo < lC\\Xe - Xl lieu 
and hence 

\\AAXe) - Ae{xl)\\ci < ^G\\Xe - X*Jci. 


Finally, for higher powers of we obtain (using (A.33)) 

||(A,)"(x,) - {A,nxt)\y < -anilx, - x:||ci , 


Tj2n 

where := C'( 2 ii)I {n€^). 


□ 


At this point we finally obtain the existence of a unique solution to (A.l) in for all hxed 
small e by applying Weissinger’s fixed point theorem f |43)1. Note that the factor 1/e in the 
estimate (A.32) provided by Proposition [A^ does not cause any trouble. Its only effect is that 
the approximating sequence {AgA{x£) converges to the fixed point slower as e gets smaller. 
Nevertheless we obtain a fixed point for every fixed (small) e: 
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Theorem A.6 (Existence and uniqueness). Consider the system (A.l) with initial data (A.2), 
satisfying ( A.3[ ), ( |A.4 ). Then for all e < eo where Eq is constrained by (A.14) and for rj given 
by (A.12) we have a unique smooth solution {u^,Ze) on [«£,«£ + v]- Moreover, and Ze as 
well as their first order derivatives are uniformly bounded in e. 


Proof. Propositions A.3 and A.5 allow the application of Weissinger’s fixed point theorem 
( |43) ) for fixed e < eo and suitable r], providing thus a unique fixed point for the operator 
on the space Xe which in turn gives a unique C^-solution Xe = {ue, Ze) on [a^, Og + r]] to the 
system (A.l) with data (A.2). Moreover, since the right hand sides of (A.l) are smooth the 


solution is smooth as well. 

The solution obtained via the hxed point argument is unique in the space Xe and thereby 
unique among all smooth solutions assuming this data by the usual argument from ODE- 
theory. 

Einally, hg, z^, and ig are bounded uniformly in e on [a^, + rj\ by the very definition of 

X,. □ 


Appendix B. Limits 


In this appendix we deal with the explicit form of the limits Ap = Zpfififi), B = 

lim£_j.o I4(/3£), and C = lim^^o as stated in Proposition 5.3 Since the actual calcula¬ 

tions are overly technical we only sketch the main points. 

Again the sign-difference between the Zj-components and Z 4 is minor and to simplify the 
notation we will use a similar convention as in Appendix and write and Z instead of 
Zpg and Zp and analogously for their derivatives. Also we will write DH instead of H^p. 
Starting with Ap we use the differential equation ( |2.17[ ) for and the uniform converge of 

to show that 


Zp£ and [4 established in Theorem 


(B.l) A = lim Zfifi,) = ^ {dH{Z^) + £2 - DH{Z^)Z^)^ + 


5.1 


zo 


Z°. 


To begin with we express Zfifie) according to (2.17) 
Z£(/3£) = 4°+ r Zfir) dr 


(B.2) 


= 4 ° + - 
^ 2 




DHdeU'^ dr - 




eZp 


-! 
2 A 


-k 


=: zO + b + IU + IIU + IVe + V^ + Vb 


U^DHdeZl 

aa^ - UfH{Zfi5, 

ilrDHZfisUeZe 


(Ta2 - UfHiZfid, 
/•A 

dr -k - 


dr 


-[ 


U^H5'fiJ,Ze 


aaf - U^H{Zfi5, 


dr -k 


f 

J Qp 


aa? - U‘fH{Zfi5, 

2 _ 

dr 


dr 


U^H5eZ, 


aa^ - UlH{Zfi5, 


where we have used that 

(B.3) \UeGe - Ue{H{Zfi5eUe) = ^U^DHdeZs - ^U^HS'.Ue - UsDHZfieUe - U^H5e 
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Proceeding term by term we have 
1 1 


\e- 


DH{Z^)U^\ = 2 ll" {dH{Z,{U-\s)))6,{s)MU7\s)) - DH{Z^)6,{s)U^) ds 


1 

< A sup 

2 weUe ^([-e,e]) 


DH{Z,{w))Ue{w) - DH{Z^)U^ 


iLi 


0 . 


where we have used that e,e]) = [ag,/?^] together with Lemma A.2, The next term, 

Me, vanishes in the limit by the uniform boundedness of the integrand, the same holds true 
for Vg. Now lllg can be treated as Ig, additionally using (5.6) to conclude 

(B.4) m, , lU^DHiZ^){Z^)^ 

We treat IV^ using f 5'(s)sds = —1 to obtain 


IV, + 


2aa^ 

1 


l7,(wmZ,(w))Z,(w) 


(B.5) 2 ([-£,£]) - U^(w)ff(Z^(w))Ss(l7e(w)) 

Finally, the limit of VI, is proportional to the limit of IV,, 

U^ff(Z^)Z° 


U^H{Z^)Z^ 




II'^£('S)'S|Ili —>0. 


vu- 


(B.6) 


aa^ 

< sup 
luef/y^ ([—£,£]) 


Ueiw)HiZ,iw))Z,iw) 


U°H{Z^)Z^ 


aa^ - U^{w)H{Ze{w))6eiUe{w)) 


aa^ 


iLi 


0 . 


By adding up the terms and using (3.10) we establish (^B.l). 


The calculations for B are relatively simple. Using equation (2.17) for I 4 we write (cf. (4.3)) 
(B.7) 


U,(/3,) = J r HiZ,irMmr))U,{rf dr ds + 0{e). 


•OLe 'JOLe 


We substitute twice, use <^c(r) dr ds = 1 and insert appropriate terms to obtain 

J rH{Z,irMiUeir))Usirfdrds-HiZ^ 

^ J a. fr J OifT 


L 


-, Ueiur^i)) J- 


HiZ,iU-\r)))6Ur)UeiUrHT))dTdl 


(B.8) 


/: 


Ue{U-\l)) 


H{Z^)5'^{t) drdl 


< 


4||p'| 


t/o 


sup 

([—£,£]) 


u,{ur\i)) J-e 

\H{Z,{w))Ue{w) - H{Z^)U^\ 

+ \H{Z^)\ sup \U,{w)-U^\ 
w>eur^([-£,£]) 


which goes to zero by the uniform convergence of Z, and U,, establishing the claimed form 
of B. 
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Finally we turn to the calculation of C which is the most demanding one. As above we express 


Ve{fie) using the geodesic equation (2.17) to obtain 


Vem = V^+ [ ^ lH{Z,{r))S',{Ue{r))U^{r)dr+ [ ^ DH{Z,{r))6e{U,ir))U,{r)Z,{r) dr 

J OLfr ^ J OL^ 

(B.9) - [ 

J a 


L 

L 


e {V,{r) + H{Z,{r))6e{Ueirm{r)) 
aa^ - H{Z,ir))6e{Ueir))Uiir) 

/?. (if7|(r)4(r) - U,ir){HiZ,ir))U,ir)5,iUeir)))) V,ir) 


dr 


aa^ - H{Z,{r))6emr))Ui{r) 

(\U^,{r)Ge{r) - Ue{r){H{Z,{r))Ue{r)6,{U,ir)))) H{Ze{r))5e{Ue{r))Ue{r) 


aa? - H{Z,{r))5e{U,{r))Ul{r) 


dr 


=: 1/0 + b + lb + llb + IVe + V, 


Note that lllg —?• 0 because the integrand is uniformly bounded. Now we rewrite 1^ substituting 

/-!/ 


s = Us{r), abbreviating w := ^(s), and using equation (2.17) for 


\^ = l- f H{Ze{w))6'^{s)Ue{w)ds = 0-l- f 6e{s){H{Ze{w))Ue{w)) ds 


(B.IO) =-- 6,is)DH{Z,iw))Z,iw)ds 


£ 

+ ^J 5MH{Z,{w)) 


\Ue{w)Ge{w) - {H{Zs)deUe)iw)j S 

aa^ - s‘^H{Z,{w))5e{s) 


ds + 0(e). 


Now the integrals on the right-hand-side of (B.IO) combine with ll^ and to give 


(B.ll) 


K(/3,) = 14°+ -!!, +Ill, + IV,+ -V,. 


Now we insert equation (2.17) for Z, into II, and follow the same pattern as before. The only 


remarkable new point is the occurrence of the regularisation-dependent term (5,(s)^s ds, 
whose prefactors cancel after a long and tedious calculation, where we repeatedly use identities 
as e.g. <^£(^) dr ds = For example we obtain for the term in (B.ll[) related to Z^ 


^ DH{Z,{r))5e{Ue{r))Ue{r)Z,{r) dr - '^DH{Z^)Z^ 


1 
2 


(DH{Z,{U-\s)))5e{s)Z,{U-\s)) - DH{Z^)5e{s)Z^) ds 


1 , 


< sup -||/5 ||li 

W&Uc ^([—£,£]) 


(\DH{Z,{w))\ 


4° - 


+ 


70 


Z' 


\DH{Z^{w)) - DH{Z^)\ ) ^0. 
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